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Abstract 


Systems with local dynamics are characterized by a finite velocity of propagation of perturba¬ 
tions, known as the Lieb-Robinson velocity. On the other hand, irreducible stochastic processes 
drive states towards some unique fixed point. However, combining both effects is mathematically 
challenging. The bounds on propagation do not depend on system-size, while the theory of mix¬ 
ing is mostly based on extensive upper-bounds. In 141 ]. a class of local Lindbladian operators on 


arbitrary lattice was constructed, for which the two effects could be combined. In this paper, we 
show that for some local dynamics, local observables are propagated inside a much smaller convex 
subset of the total Banach space. This allows us to show localization for such dynamics. 
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I. INTRODUCTION 


In 1958, Anderson [8] showed that random impurity potentials produce the localization 
of quantum wavefunctions. This lead to a better understanding of transitions between 
insulators and conductors. In quantum spin-lattice, it can be shown that due to (quasi)- 
local interactions perturbations propagate at a finite velocity. This was proven by Lieb and 


ifl. 


Robinson [If. Recently, there has been an increasing interest in studying local dissipative 
dynamics. One of the major drive behind, is the possibility of using dissipation as a tool for 
engineering states. Similarly to classical Monte-Carlo algorithm, the goal of course is to drive 
any state as fast as possible towards a desired state. When introducing local dissipation to 
the system, the propagation of the perturbation competes with the global relaxation of the 

a, A, 


system. In 


a study of such competition was first attempted. In 


14], it was shown 


that this competition could not only effectively reduce the velocity of propagation but also 
lead to some increasing deceleration. All the results so far were studied for systems with 
the uniform state as unique fixed point. For some development of quantum technologies, 
localization of excitations on pure states such as graph state is of general interests js]. 


A. Local Quantum Stochastic Dynamics 

Consider a D -dimensional lattice Z D with a metric. At each point x E Z D of the lattice, 
define a d-dimensional Hilbert space T~L X and for each finite set A C Z D denote the product 
space, 

U A = (g)U x 

Denote the algebra of all matrices acting on "Ha, he. the algebra of local observables of A, 

A a , 

A a = ®B(H x ) 

xgA 

If Ai C A 2 , the algebra A Al can be identified with the algebra A Al <g> 1 ai\a 2 and therefore 
A\i C A A2 . Define the support of a local observable A E A a as the minimal set X C A for 
which A = A' ® l A \x- 

The dynamics of the systems is generated by some local Lindbladian, 

dtA(t) — L[A(t)\, L[.] — I z [] — B} z []R z + Qz[] + [-]Qz> ^ R\Rz + Qz + Q\ = Qz £ Az 

z z z 
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where the last condition insures that the dynamics preserves unity. 


B. From local Response to Lieb-Robinson Velocity 

By Frobenius-Perron, the adjoint of the dynamics has some positive fixed point p eq , 

L O d [Peq] = 0, Peg > 0, Tr p eq = 1 

Let us now perturb the dynamics at some local subset Z, be it time-independently or not, 

dtA(t) = L 0 [A(t)\ + (f)z{t)[A{t)\ 

We wish then to compare the expectation value of the undisturbed system (A(t)) = (A) 
with the new one, (A(t))$. This is readily derived, 

(A(t))^ ~{A) = f Tr (p eq 4> z [A(s)]) 

Jo 

Hence the response of the systems, is determined by bounding the term inside the integral. 
Since Tr(p eq ) = 1, this becomes, 

Tr (p eq( j) Z [A(t)]) < \\(j) Z [A(s)]\\ 

Such expression was studied for the first time by Lieb and Robinson and by reducing the 
problem to a Laplacian inequality it was shown that, 

||fe[/lW]||<Cexp(^ i ' t ~^ j4) ^ (1) 

where C , v and £ are finite constants. 


C. Killing of Transport 

In [l^], it was shown that under some symmetry a class of Lindbladians could be con¬ 
structed so that the bound takes the form, 


Uz[A(a)]\\ < C exp 


f* exp (-ps)v ds — d(Z, A) 

e 


( 2 ) 


As pointed out, if one could show, 

exp (tL c ) [X] - lim^oo exp ( tL c ) [X]\\ 


max ■ 
x 


\\X\\ 


< C"exp(—A t) 


(3) 
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where C c is the full generator of the dynamics minus some local interaction terms, and C 
is a finite constant, the bound can be shown to be of the form, 

II^Wlll < e"exp ( «*-»-^M) ) (4) 

This becomes particularly interesting when A ~ v, as this imply that A(t) is localized inside 
a fixed volume for all t. Proving the case for which an inequality such as fl3]) is satisfied is 
hard. We present here an alternative. We show a result similar to d3J), by relaxing the set 
of (global) observables to a smaller convex subset. Under the additional property that the 
dynamics preserve the convex structure, a bound such as (j3J) can then be proven, and even 
more so localization can be implied. 

II. LOCALIZATION 


Consider some local basis {cr^} of the local Hilbert space C^. For the lattice aC^, 
define the convex hull of the set constructed from the local basis multiplied by some phase 
C = Conv <S>j cr^\(p G [0, 27r]}j. We say that this set is normalizable with local observ¬ 
ables, if VA e Aa a , there is some hnite constant 0 < Ca < oo, so that, A/Ca G C. 

We say that a local dynamics preserves C if, 


Vt, X e C, 3 Y(t) EC, 0 > r(t) < 1, exp (tL)[X] = r(t)Y(t) 


(5) 


As, we see further on, it is then sufficient to restrict ((3j) to the inequality, 

exp (tL c ) [X] - lirn^ooexp ( tL c ) [X]\\ 


max 

XGC,||X||=1 


ll^ll 


< exp (—At) 


( 6 ) 


where C c denotes the global generators minus any number of local interactions Iz- As usual 
in order to take into account long range interaction, we introduce the usual formalism of 
reproducing functions. We assume there exists a non-increasing function F : [0, oo) —y (0, oo) 
so that F is uniformly integrable over the lattice, 


Ill'll := sup F ( d ( x ,v)) 


< oo 


x£Z c 


y&z c 


and there is a constant C, 




F {d(x, y))F(d(y , z) ^ d ( x , y )+d(y,z)-d(x, z )) 
F(d(x, z)) 

zGZ a 
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so that C p | < oo. 

Assume there exists a /x > 0, so that, 


\\ L \\n= sup 


1 Z\\cb 


?ud{x,y) 


x,yez d z ^y F ( d ( x ,y)) 


< oo 


Finally, let us assume there exists some finite constant 0 < C p hi < oo so that VZ, WX G C, 

L z [X}/C x eC, C x < CfdX 
We say that the system is frustration free if, 


VZ, Cz ° lim exp(f£) = 0 

> co 

In order to derive localization, the dynamics needs to be split in two other local one. 
The first one is responsible for the spreading the perturbation, while the second is the 
main dissipative part which erases the perturbation. Clearly, for localization to be possible, 
perturbation needs to be driven towards some zero-correlated subspace. Hence, we assume 
that the local parts of C\ = Ylz Iz mutually commute. For the sake of the fluidity of the 
argument, the interaction terms of both generators are grouped inside some finite subset Z. 

l = lq + Li = y^(/|+ if) 

z 

We can now state, the theorem. 


Result 1. Let L = L 0 + Li = + -Tz) frustration free and irreducible local Lindblad 

generator satisfying (EP and (EJj, and [If, If,) = 0. Then VA G A\ a , B G A\ b with D A B = 


||[R,exp(fL)[A]]|| < 2 


Ill’ll 


C„ 


\\B\\ cb min{A A , A B } (exp [vt - At] - 1) exp [-/id(A, B)) (7) 


with v = and, 


max H^P(iU)[A1-lim, expftLQI.Ylll < 
Xec,\\x\\=i X 


( 8 ) 


Hence, we see that the velocity and the decay rate are decoupled from one another. By 
varying the way between the diffusive part £° an the dissipative part £ 1 , it is then possible 
to get A > v, from which localization follows. 

As the proof of the result is extremely similar to the proof presented in the ” further 


discussion” of 


14], we leave the proof in the appendix. 
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A. Classical Stochastic Processes 


It is straightforward to embed a classical stochastic dynamics into a quantum one. For 
example in the case of a local stochastic generator L = 'Yf z Tz — Iz, where Iz is the identity 
operator and T z is local stochastic. Hence, this can be embedded into a Lindbladian using 
the local standard basis { e(Z)j }, 

41 = E E \ATizhZ) <Z)AWZ)] - {.}z (9) 

^ 3 

In the case of local Hilbert spaces c|, e(Z)ij = \i)(j\ . Additionally, we can see a clear de¬ 
coupling between the vector spaces span{<8crj|cr,' G {l, a z }} and span{<g)<Tj|cr,- G {1, a x , cr y }}. 
Hence for simplicity, let us restrict to dynamics on lattices <g)j € aC 2 , and dynamics (0 re¬ 
stricted to observables in span{®crj|crj G { 1 , cr^ }} 

Define now the local unitary Uj = ( a z + a x ) /\/2. It is then also understood that Uz = 
®j e zUj. In the following definition the matrix \A\ iy is understood as the entry-wise absolute 
value. The following proposition gives us a local condition that implies property (J5]). 

Proposition 2. If the matrix \UlTllIz\ij is sub-stochastic then the dynamics (G|) preserves 
the convex set span{e*^ <g) cr^|cr^ G {1, cr 2 }, cj) G [0, 2n]} as defined in $5\). 

Proof. Since, 

exp(A + B) = lim ((/ + A/n)(I + B/n)) n 

n—>■ oo 

and due to the dehnition of ([9]), we can decompose each interation term as (1 — e)/ + eT. 
Hence, it is sufficient to show that, \/Tr(e(Z)jT z ) e(Z)j[.]e(Z)j satisfy the claim. It can 
be seen that the claim is equivalent to showing that Vr G {^crj\aj G {l,cq}} , 

£ \{o\T z \t)\ <2' z ' 

IcTj e{l,cr z }} 

which is indeed equivalent to \U z T z Uz\ij being sub-stochastic. □ 

As an illustrative example, we can consider the dynamics, 

41 = £“[.]+ A [.]-[.] 
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Clearly we can see, 


max || exp ( tL c ) [X] — lim exp ( tL c ) [X]|| < lim ||(1 — t\/n)[X\ + t\/2n'S~' TrJX] — P[X]]|| n 

C,||-Y|| = l t—>oo n—too J 

j 

< exp (—At) 


with P[X\ = (i) |A| . 


B. Engineering dynamics of Graph States 

As a second and final illustration, we show that some Lindbladians known to generate 
graph states satisfy the properties (JHJ) and (J5]) for C = spanje*^ <g) <7j\<7j G {l, a x , a y , a z }, (j) G 
[0,27r]}. Graph states are defined as the unique eigenstate of the set of commuting ob¬ 
servables Uk = cr T y. ®j<zz cr z .j, with Z some neighborhood k. As shown in Q, defining 
c k = + Uk)<J x ,a with a G {x,y}, the graph state can be prepared by each local Lindbla- 

dian, 

£“[■] = J2 c k ] ^ c< k - \{ c k tc k,-}, « e {x,y,z} 

k 

Notice that P k = + £4) is a projector. Let us therefore take £ G {<Z<jj\(jj G 

{l, a x , a y , <J Z }}- Let us fix a = x as the other case are proven similarly. We need then 
to distinguish between the cases a a ^£ = ±£<70,^ and £4£ = ±££4 and calculate for each, 

1. Ua k'Z = £<J Qi fc, UkZi = ££4 implies that ££[£] = 0 and so exp(LCj£)[£] = £. 

2. cr Qt fc£ = £cr Q!jfc , U k T, = —T,Uk gives ££[£] = —£ and so, exp (££*:)[£] = e _At £ 

3. cr a)fe £ = —£cr a fe , £/ fc £ = TJJ k yields ££[£] = -£ - £cr Q)fc C4cr Qifc and so, exp(f£^)[£] = 
£-Pfc + e _2At £(l — Pk) 

4. cr aj fc£ = — £cr a) fc, f/fc£ = —££4 also gives ££[£] = —£ and so, exp(££j£)[£] = e _At £ 

Since cr a ,kUkCra,k G C, all cases therefore imply ([5]). Finally in order to prove (JHJ) , since per 
construction the Lindbladian has a dimensional one kernel, case 1 cannot be satisfied for all 
k. The claim also follows if either case 2 or 4 are fulfilled for a single k. Lienee, the only 
remaining point to look at, is when case 3 is satisfied at least once. So for all k, [£, P k \ = 0, 
and we can take mutual eigenvectors |a), £|a) = s a |a), sa = ±1 and P k \oi) = p a ,k\oi), 
Pa,k G {0,1}. Clearly, since the graph state is per definition the unit vector which has 
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eigenvalue 1 for all Pk, it has to be an eigenvector of E too. Hence for all |a) orthogonal to 
the graph state, we have, 

| (a| exp(f£“)[£] - lim exp(iC“)[£]|a)| < e~ 2t 

t—yoo 

Hence all case combined, we see that the decaying condition (jHJ) is satisfied. 

III. CONCLUSION 

We studied a class of local Lindbladians which preserves a particular local convex struc¬ 
ture. We show that in the irreducible case, a competition can emerge between the diffusive 
and strongly mixing part of the generator. We see that the diffusion property can be decou¬ 
pled from the mixing. Localization is then implied for a finite mixing rate of the same order 
as the Lieb Robinson velocity. At the end, we showed various construction from classical 
stochastic processes, to more complex quantum dynamics with graph states as fixed point. 


Appendix A: Proof of Result 1 


Proof. The proof is very similar to the proofs presented in [REF DESCAMPS] Define a path 
of n steps, V n , between A and B as the set of the support of the local interactions I\ j} 

V n = {Aj|A# n Ai 7 ^ 0,Ai n A 2 7 ^ 0,..., A n D A a 7 ^ 0 , all other intersections are empty} 


Denote V* as the subset of V n containing the first k steps of the path, where V a n — 0. Define 
the generator part C° containing local interaction which do not intersect with the path, 


L k = L ~ E '' 

znv&® 


In Lemma 1 we show, 


with, 

JVn = 


71=1 


t pt-S 1 /**-£”=! Sj 


||[fl,exp((LM]||<^||L||R 3 upJ I .„)^ ]T n E E 

' Vn ' V n x 0 eA B Aj&VnXjGAj x n+ 1 eA A 



0 ^0 


ds [B ,.] ^exp 

3 = 1 


s ]L v j-i 

I n 




A,-/ \\ 1 A j \\cb 


ex P[(^- J2 s ^ lc ^a 

3 = 1 


(Al) 








with the ordered product IIj=i Cj = C\C 2 ---C n and the supremum is taken over all paths of 
length n. 

Dehne the projector, 

P = lim exp(t£ 1 ) 

t—>oo 

The frustration free and irreducible condition allow us to insert the projectors, 


Jv„ - 


t rt.-S! rt~127= l s j 



0 Jo 


ds [B,.] (exp 


S jL-p:i-i 

1 n 


(i-^vunju 


3 = 1 


ex P[0 - J2 s j) L r«\ (! - p ) A 

3 =1 

Since, [L%, kl Pf, k ] = 0, \/V n ,k and using the assumption (JSJ), it follows, 

'n 'n 

ClWLXt" 


Jv,_ < exp(-At)- 


n\ 


From which follows, 


||[Ar t ( J B)]|| < C'\\A\\ ||f?||exp 


-At + vt — d(A, B ) 

I 


□ 


Lemma 3. For A and Iz with Z D = 0, 

||/Iexp((L)Al|l<X;i|I“l|Asup* E 11 V v FXx^XjOe-^-^ 

n =1 ' Vrl ' Pn XoGAa Aj£Vn XjGAj Xn+lGZ 

Proof. The idea of the proof, is the gradually take out, each interaction term of the diffusive 
generator £° forming so a path from Z to A. Define the generator 

r° _ V T° 

^nz ~ 1 A 1 

AiDZ^0 

We can then expand, 


I z exp(tL)A = ^2 II^aJI cb f IzexplsoL^JWlAAcbexpft - s 0 ){L - L° nZ )\A 

Ain z^H 

We then continue this procedure for each term Ja x exp[(t— sq)(L—L ff z )]A for which AiflA^ = 
0. Let us say that after k — 1 steps we have A^ D A a ^ 0, from our construction we see that 
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we have built a path of k steps from Z to A for which 3v k {z,A) can be bounded by equation 
(ED. For each n, we can then use the definition of HL^, 

fsup J Vn ) Yi n ^h j aJ c6 - ( su p^) sup n s F ( d ( x j-h x i)) e 

^ Vn ' Tn AjGVn Aj \v„ J Vn xoGAa AjGVnXjGAj x n+1 GZ 

□ 
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